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String breaking by dynamical quarks in (2+l)-d lattice QCD is demonstrated in this project, by measuring 
the static potential and the local color-electric field strength between a heavy quark and antiquark pair at large 
separations. Simulations are done for unquenched SU(2) color with two flavors of staggered quarks. An improved 
gluon action is used which allows simulations to be done on coarse lattices, providing an extremely efficient means 
to access the quark separations and propagation times at which string breaking occurs. The static quark potential 
is extracted using only Wilson loop operators and hence no valence quarks are present in the trial states. Results 
give unambiguous evidence for string breaking as the static quark potential completely saturates at twice the 
heavy-light meson mass at large separations. It is also shown that the local color-electric field strength between 
the quark pair tends toward vacuum values at large separations. Implications of these results for unquenched 
simulations of QCD in 4-d are drawn. 



In unquenched simulations, we expect the lin- 
ear rising static quark potential to saturate at 
large quark separations R. This is the phe- 
nomenon of hadronic string breaking. However, 
despite extensive simulations of full QCD by sev- 
eral large scale collaborations (l), a convincing 
demonstration of string breaking remains contro- 
versial. Traditionally, Wilson loops have been 
used to extract the static quark potential, though 
recently it has been suggested that this operator 
is not suitable because it has a small overlap with 
the two- meson ground state at large R . 

On the other hand, the Wilson loop produces a 
heavy quark-antiquark trial state without explicit 
light valence quarks. This trial state is of great 
physical interest because of the analogy with the 
process of hadronization, where light quarks are 
not present in the initial state, but rather are ma- 
terialized from the vacuum. 

It was proposed in Ref. ^ that string break- 
ing can in fact be observed in Wilson loops by 
working on coarse lattices using improved actions. 
On coarse lattices the computational effort can go 
towards generating much higher statistics, which 
allows much longer length and time scales to be 
accessed. This approach was applied in Ref. [D 
to lattice QCD in (2+1) dimensions, and later to 



4-d QCD in Ref. [Q. It now appears that the 
overlap of the Wilson loop with the ground state 
is large enough, in the relevant range of R, to 
allow string breaking to be resolved. 

The present work is a follow-up on the earlier 
study of Ref. [||. What is new here is a much 
more extensive study on larger lattices and with 
far higher statistics. The results give a much more 
convincing demonstration of an asymptote in the 
static quark potential at large quark separations. 
The demonstration of string breaking from mea- 
surements of the local color-electric field strength 
in unquenched QCD is entirely new to this work. 

To reduce computational cost, we work in 3-d 
unquenched QCD with SU(2) color. QCD3 has 
been shown to share most of the fundamental 
properties of the realistic 4-d QCD Q and thus 
provides excellent insights to QCD. The coupling 
(7q in QCD3 has dimension of mass and so ex- 
plicitly sets the mass scale for the theory, i.e., 
m cx (7g for any mass quantity m. In order to 
have some intuition for this mass scale, we ex- 
press this scale in terms of the length scale in 4-d 
— "fermi" (fm). Motivated by the fact that both 
3-d and 4-d QCD are linearly confining, we match 
the quenched continuum extrapolation of the 3- 
d string tension, ^/a/g^ = 0.3353(18) (|, with 
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the physical value of the string tension in four di- 
mensions, = OAAGeV. This in effect sets the 
value of go for our 3-d theory. A similar value 
for the coupling is obtained from other physical 
quantities, such as equating the unquenched p 
meson mass in this theory with the physical value 
in 4-d QCD. 

The tree-level 0(a^)-accurate improved gluon 
action in 3-d has the usual form 
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where P^^ is the plaquette and R^i, is the 2x1 
rectangle. The bare lattice anisotropy is entered 
through = as/at and ^ij = at /as = 

1,2). We use the staggered quark action in 3-d, 
which has the same form as the 4-d action and 
describes two flavors of four- component spinors: 



Sks = E Ct^'np.{x)x{x) [U^,{x)x{x + fi) 

X.fJ, 

-Ulix - fi)x{x - A)] + 2amo ^ xix)x{x), (2) 

X 

where ^3 = ag/at and ^1^2 = 1- 

For comparison, most of the simulation param- 
eters employed in the present project are taken 
from Ref . |^ . Simulations were done on a 22^ x 28 
lattice with /3 = 3.0. Using the scale setting 
procedure discussed earlier, the lattice spacing is 
identified to be ~ 0.2 /to in physical units. 
The input bare quark mass in lattice units is 
scaled according to TOo/fifg = 0.10, with m^/mp 
found to be ~ 0.75. The input bare anisotropy 
is at /as — 1/2. The configurations were gener- 
ated using the HMD algorithm (the ^-algorithm) . 
Overall 30,000 measurements were taken, corre- 
sponding to a total of 300,000 trajectories. 

Only (fuzzy) Wilson loops were used to ex- 
tract the static quark potential. The final results 
are presented in Fig. for propagation times 
T/as = 2,6,12 (or T ~ 0.2, 0.6, 1.2/m in phys- 
ical units). The unquenched heavy- light meson 
mass was also computed and is shown as the dot- 
ted lines in Fig. |l|. The results here give clear 
indication of string breaking as the unquenched 
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Figure 1. The static quark potential V{R,T). 



potential is substantially flattened at large R and 
approaches the expected asymptotic value of the 
heavy-light meson mass. 

To emphasize the flattening of the potential 
at large values of R, we fitted the data points 
within R/as = 8.0 — 10.0 with a linear function 
V{R) — aR + b and plotted the slope a (string 
tension) against the propagation time T in Fig. 
||. Results here unambiguously show that the 
slope tends toward zero at sufficiently large T, 
T~ 1.2/to. 
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Figure 2. String tension ala as a function of T. 
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Figure 3. Flux tube profiles for R/qs — 5,7. 
The open squares and the solid circles are the 
quenched and unquenched results respectively. 



The lattice observable needed to measure the 
field strength is given by the following correlator 
of plaquette P^i, with a Wilson loop W{R, T) [|| 

fR.,Ti^. <W{R,T) (P^.(f)-P^.(foo)) > 

<W^(i?,T)> ' 

Here, x is measured relative to the center of 
the Wilson loop and x^o is usually taken as the 
site half the lattice away on a finite size lattice. 
We measured the /13 component, which is cor- 
responding to the energy density of the color- 
electric field in the direction parallel to R in the 
continuum limit, i.e., /13 = The flux tube 
profiles for quark separations R/ag = 5,7 are 
given in Fig. ^ where f "(a;) is plotted both along 
{x\\ ) and perpendicular {x±) to R. The data is ob- 
tained with T/os = 7. Results here illustrate the 
formation of flux tubes between the color charges. 
Notice that the thickness of the tube is ~ 1.6 fm 
and is independent of the quark separation R. 

Fig. ^ also demonstrates the screening of 
the sources by the dynamical quarks: the un- 
quenched field strength is significantly lower than 
the quenched results inside the flux tube. One 
would also expect the unquenched field strength 
to drop substantially while the quenched results 
stay fairly constant near the string breaking dis- 
tance. This idea is illustrated in Fig. where the 



Figure 4. Unquenched i?"(0) for different R. 



field strength at the center of the tube is plotted 
against the quark separation for the unquenched 
data. One can observe that at sufficiently large T, 
the field strength decreases gradually when R in- 
creases, and eventually vanishes beyond the string 
breaking distance of R/ug ~ 8.0 — 10.0. 

To summarize, we were able to observe string 
breaking by dynamical quarks in QCD3 . The use 
of improved actions on coarse lattices provided 
a crucial advantage in accessing the quark sep- 
arations and propagation times at which string 
breaking occurs. The results confirm that string 
breaking appears in Wilson loop correlators only 
at sufficiently large T, of about 1/m. These have 
clear implications for simulations of full QCD in 
4-d as this scale is computationally accessible par- 
ticularly on modestly coarse lattices. 
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